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Abstract. This letter examines the consequences of a recently proposed modification 
of the postulate of equal a priori probability in quantum statistical mechanics. This 
modification, called the quantum microcanonical postulate (QMP), asserts that for a 
system in microcanonical equilibrium all pure quantum states having the same energy 
expectation value are realised with equal probability. A simple model of a quantum 
system that obeys the QMP and that has a nondegenerate spectrum with equally 
J> I spaced energy eigenvalues is studied. This model admits a closed-form expression for 

' the density of states in terms of the energy eigenvalues. It is shown that in the limit as 

. the number of energy levels approaches infinity, the expression for the density of states 

' converges to a J function centred at the intermediate value (i?max + £'min)/2 of the 

■ energy. Determining this limit requires an elaborate asymptotic study of an infinite 
, sum whose terms alternate in sign. 

:^ 
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1. Introduction 

This letter investigates a generalization of the usual definition of a quantum system 
in microcanonical equilibrium. If the Hamiltonian H that describes a system has 
a nondegenerate spectrum, then according to the standard definition of quantum 
microcanonical equilibrium the system must be in one of the eigenstates of H. This 
requirement is known as the postulate of equal a priori probabilities [Ij. We emphasize 
that according to the definition of microcanonical equilibrium in Ref. [Tj the state of 
such a system cannot be a linear combination of eigenstates of H. However, if H has 
a degenerate spectrum, then the density matrix that describes a system of energy E in 
microcanonical equilibrium contains all states \E,k) of the degenerate energy E with 
equal weight: ^ ^^=1 \E, k){E, k\, where n is the number of states having energy E. 

Because the standard definition of quantum microcanonical equilibrium only allows 
the system to have energies that are eigenvalues of H, an alternative, less restrictive 
definition has recently been introduced j2]. By this latter definition, called the quantum 
microcanonical postulate (QMP), a state of a system in microcanonical equilibrium can 
have an energy that is not an eigenvalue of H. The discussion in Ref. j2] of quantum 
systems obeying the QMP is qualitative. Here, we give a quantitative analysis of a 
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quantum system described by a Hamiltonian having a nondegenerate, equally spaced 
spectrum. Assuming that the system is in microcanonical equilibrium and that it obeys 
the QMP, we study the behaviour of the density of states fi{E) as the number of energy 
levels becomes large. 

This letter is organised as follows: In Sec. El we review the representation for the 
density of states f^{E) in terms of the energy eigenvalues as outlined in Ref. We then 
define the model investigated in this letter in which the energy spectrum is taken to be 
nondegenerate and to grow linearly: Ek oc k. In the next two sections we investigate 
the behaviour of /i(-E) as the number of states of H becomes infinite. In Sec. |3]we show 
that fi{E) integrates to unity. Section IH presents an asymptotic study of yu(-E) as the 
number of states becomes infinite. On the basis of the analysis given in Sees. ElandlU 
we conclude that fi{E) approaches 6[E — (-Emax + -E'mm)/2]. 



2. Definition of the model 

Let us review briefly the general mathematical framework proposed in Ref. ^ for 
describing the density matrix of a mixed state of a quantum system in microcanonical 
equilibrium. Consider a quantum system defined on an {n + l)-dimensional Hilbert 
space H. Let Z°' (a = 0, 1, 2, ... , n) be a typical element of H and let denote the 
Hamiltonian with eigenvalues Ei {i = 0,1,2, ... ,n). Then, the expectation value of the 
Hamiltonian in the state Z" is (H) = Z^W^Z^ / Z.^Z'^ . Assume that in microcanonical 
equilibrium all states satisfying the condition (if) = E are realised with equal 
probability. Then, the corresponding unnormalised density of states is 

n{E) = - J (T+^Z (T+^Z S{ZaZ'' - 1) 5 f — ^ j • (1) 

The constraint — 1) in (P) arises because one is only interested in the 

unit normalised states, and the factor of vr reflects the additional redundant overall 
phase of the state. It is convenient to use the standard integral representation 
S{x) = ^ J^dAe"''^^' for each of the 6 functions appearing in ((T)). The Hilbert-space 
integration then becomes Gaussian in the Z variables leaving the expression 

n^E) = i(-i7r)"+i r du dx A 

^ ^ TT^ ^ y-oo27ry_^ 27r fi ^ + '^^^ 

Assuming that the energy spectrum is nondegenerate, one can perform the A 
integration to obtain 

The remaining u integration can now be performed explicitly to give 
n{E) = {-nrY,6'^-\E,-E) n Y-^^ 

A:=0 l=0,^k '■ ^ 
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where 6^ "^(x) denotes the nth integral of the 6 function: 




(5) 



The density of states Q{E) as defined in is normahsed by dividing it by the 
volume of the subspace of 7i spanned by states having unit length: ZaZ"' = 1. This 
gives the normalised microcanonical state density function fi{E). The volume is given 
by 7r"/n! (see, for example, Ref. |SI). Thus, n{E) = n! vr"" gives the density of 



In this letter we propose a particular QMP model in which the energy spectrum 
rises linearly and is given by Ek = k. Our objective is to study the behaviour of IJ^{E) 
as the number of energy levels becomes infinite. With this linear choice of spectrum the 
normalised density of states becomes 



where the notation [E] indicates the largest integer less than or equal to E. 

It is now convenient to rescale the energy spectrum so that the range of the energy 
lies in the interval [0, 1] for each n. Upon rescaling, (jH)) transforms to 



where E G [0, 1] for all n. To derive this result we have used the fact that the sum in 
dni) vanishes when the summation range is taken from k = to k = n. 

In Fig. HI we plot the density of states fi{E) in ((Zj) for n = 3, 6, and 9. This graph 
suggests that fi{E) converges to a 5 function centred at = 1/2 as n, the number of 
energy levels, increases. We show analytically that the density of states ^{E) associated 
with a quantum system having the spectrum Ek oc k does indeed approach S{E — 1/2) 
in the limit n — > cxd. Our analysis is of interest because it involves an asymptotic study 
of an infinite sum whose terms alternate in sign. To overcome the difficulties associated 
with this alternating series, we convert the series to a double contour integration whose 
asymptotic behaviour is obtained using the method of steepest descent. This work also 
provides a new limit identity for the Dirac 6 function. 

3. Analysis of the model 

To verify that fi{E) approaches 6{E — 1/2) as the number of states n approaches infinity, 
we must establish two properties of fi{E). First, we must show that dE fi{E) = 1. 
Second, we must show that the limiting value of yu(-E) is zero except at E = 1/2 where 
it tends to infinity. In this section we show that the normalisation condition satisfied 
by fi{E) in is valid. Let us define J by / = dE n{E). The summation in must 
be evaluated piecewise because of the dependence of the summation range on E. Thus, 



states that satisfies the normalisation condition dE fj,{E) = 1. 




(6) 




(7) 
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Figure 1. The density of states /in(£') associated with a quantum system having a 
hnear energy spectrum iS^ = k, where the range of the energy is suitably rescaled so 
that E hes in the range [0, 1] for all n. Plots of iJLn{E) are given for 4-, 7-, and 10-state 
(n — 3, 6, and 9) systems. Observe that as the number of energy levels increases, 
the distribution becomes more peaked at the centre E = 1/2, suggesting that as the 
number of energy levels approaches infinity, the distribution approaches 5{E — 1/2). 
The analysis in Sees. |2| and 01 verifies that this is indeed the case. 



it is convenient to decompose the integration range of / into n intervals and to write 



\nE] 



■lf{k-nEY-^ 
k\{n-k)\ ' 



(8) 



j=i fc=o 

To perform the integration over E we rewrite the summation in the integrand so 
that it is independent of E. Given that j ranges from 1 to n and that k ranges from 
to n — 1 with k < j — 1, we have 



dE- 



j=i k=o 



-lf{k-nE) 
A;!(n- A;)! 



n—1 



^ kUn - k)\ ni-n) 
j=l k=0 ^ ' ^ ' 



(9) 



E={j-\)ln 

i-1 



We now interchange the order of summation according to X]j=i X]l;=o ~ Sfc=o Sj=/c+i- 

71-1 < ^\k " 

" (10) 



k=0 ^ ' j=A:+l 

Performing the sum over j, we obtain 



E [(^' - - 0' - ^ - 1)1 



k\{n-k)\ 



{n - ky 



{k-nY. 



A;=0 ^ ' k=0 

Observe that the summation in is the nth discrete difference of fc". Recall 
that for the polynomial f{k) = /c" + lower powers, the first discrete difference is 
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Vf{k) = f{k) — f{k — 1) = nk^~^ + lower powers. The second discrete difference 
is V'^f{k) = f{k) — 2f{k — 1) + f{k — 2) = n{n — l)/c"~^ + lower powers, and so on. The 
nth discrete difference is especially simple because there are no remaining lower powers: 
T>^f{k) = n\. This observation allows us to evaluate the sum in : 

n 

E Ck) = ^^^^ 

We have thus verified the normalization condition 1 = 1. 



4. Asymptotic behaviour of (|7|) for large 



n 



We now examine the behaviour of /i(-E) in the limit as n — oo. We have already shown 
in Sec. El that the integral of f^{E) is unity. To show that fi{E) approaches a delta 
function as n ^ oo we must establish that yu(-E) becomes singular at the central value 
E = ^(-Emax + -E'min) and that it vanishes at all other points in this limit. Note that by 
the scaling used in the central value is at = 1/2. 

The representation of /i(-E) in ((7j) for finite values of n is symmetric about the point 
E = 1/2. To verify this symmetry we make the transformation E ^ 1 — E and replace 
the summation variable khj n — k. Thus, we need only study the behaviour of fi{E) for 
E = 1/a, where a > 2. Without loss of generality, we set n = a J, where J is a large 
integer, and let iOj{a) be the value of fi{E) ai E = 1/a: 

It is straightforward to find the behaviour of ujj{a) for large J when a > e. Using 
Stirling's formula for the asymptotic behaviour of the factorial function, we observe 
that each term in the sum in ()13|) is exponentially small; that is, it has the form e'^"^ 
(J — > oo), where A is a positive constant. The number of terms in the sum grows 
linearly with J. Thus, the sum vanishes as J — oo. 

However, when 2 < a < e, the terms in the sum (fT^ are exponentially large. In 
this case, the factor of (—1)'^ in the summand gives rise to a deep global cancellation 
among all the terms in the sum. When a > 2, this cancellation causes ujj{a) to vanish 
exponentially for large J. The case a = 2 is special because the sum does not vanish 
exponentially. We have performed the sum on the right side of fll3|) numerically for 
large values of J when a = 2 using Richardson extrapolation [4 . We find that 

cjj(2) ~ (1.9544100476... (J ^ oo). (14) 

Establishing these asymptotic results analytically is difficult. Laplace's method for sums 
cannot be used to evaluate ujj{a) because Laplace's method involves local analysis and 
this method is inadequate when terms in the sum alternate in sign. 

To overcome this difficulty we convert the sum in (fT^ to a double complex contour 
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integral. We begin by substituting k = J ~ p: 

^ T{J -p+l)T[{a -l)J + p+l] 

We then use the identity = f^dteH~^ to represent the T functions in (fT3j) . The 
contour C is infinite and encloses the negative real-t axis; C can be taken to be a circle 
around the origin when a is an integer. Rescaling the integration variable t gives 



I drdsr-("-i)^-is-^-^y-f-V+^y, (16) 
(27ri)2 JcJc^ Zr[P^ r J 



where we have interchanged orders of integration and summation. 
Integrating by parts with respect to r, we obtain 

■'^ ' ^ ' (27ri)2/c/^, r s ^ ^ 

We also integrate (fTB|) by parts with respect to s: 

a^J f f dr ds 

ujjia 



'CJC ■ 



We then evaluate the finite geometric sums in (|T7jl and (|T8|l using the identity ^p^i = 
(a"^"*"^ — a)/{a — 1), where a = —e^e^s/r. The representation for ujj{a) simplifies when 
we combine the right side of (fTTj) multiplied by (a — l)/a and the right side of (fTHj) 
multiplied by 1/a, and then replace s by —s: 



(2vri) 



drds (s-r)r-"-^-^e" . (19) 



The term proportional to e"'*^''"*'' in the integrand of (fTI?|l is analytic in s along the 
real-s axis for s < 0. Hence, by shrinking the contour to a small circle about the origin 
in the complex-s plane, we find that the integrand does not contribute to the asymptotic 
behaviour of for large J. We have thus reduced the expression for LJj{a) to 

Also, because the integral fdss^"^ vanishes for integer J > 1, we may simplify ()2()j) 
further by adding s~"' to the integrand of the s integral: 

Note that the integrand of pTjl is singular if 

se^ - re' = 0, (22) 



as long as e''— does not vanish. Clearly, (j22|l is satisfied when r = s, but the numerator 
of the integrand in (pT|l also vanishes when r = s. Thus, it may appear at first that there 
is no singularity at r = s. However, for the special point r = s = 1 the denominator 
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has a higher-order zero than the numerator and thus the integrand is singular there. 
To find the asymptotic behaviour of ()21|) we must perform a steepest-descent analysis. 
However, if we look for a saddle-point of the double integral, we find that it is located 
near the singular point r = 1 and s = 1, which complicates the asymptotic analysis 
enormously. Instead, we will evaluate the s integral in closed but implicit form and 
evaluate the remaining single integral in r asymptotically. 

It is remarkable that the transcendental equation (j^^ has other solutions for which 
r ^ s. These solutions cannot be expressed in closed form. However, we have discovered 
an explicit parametric solution to ()22|) for which r ^ s: 

r = X e~^/ sinh A and s = A e^/ sinh A, (23) 

where A is any complex number 

To evaluate the s integral in closed form we must take the asymmetric solutions 
(I23|l into account. We treat the C contour as a circle about the origin in the complex-s 
plane, but rather than considering the singularities inside this circle, we include instead 
the contributions of the singularities outside this circle because the integrand vanishes 
at |s| = oo in all directions. We now solve ()22|) for s as a function of r and denote the 
solution as s = S{r). We then use residue calculus to evaluate the integral (PT|) at the 
simple pole located at s = S{r). The result is 



where we have simplified the integrand by using the algebraic relation in 

To prepare for the asymptotic evaluation of the integral in (j24j) we rewrite it in 
standard Laplace form in terms of the parametric variable A in 



u;j{a)= d>dXg{X)e-'^^'\ (25) 



where /(A) = log (Ae^/ sinh A) + (« — 1) log (Ae ^/ sinh A) and g{X) = a J sinh A(l — A — 
A/tanhA)/[i7r(sinhA - Ae^)]. 

Following standard steepest-descent techniques we identify the saddle point as 
the solution Aq to /'(A) = 0, where /'(A) = 2 — a + a{l/X — 1/ tanhA). It is easy 
to verify that /(Aq) > when a > 2. This implies that ujj{a) vanishes exponentially 
rapidly like ujj{a) ~ e~^^^°^'^ as J — > oo for a > 2. However, when a = 2, Aq = and 
/(Aq) = 0. In this case ujj{a) behaves algebraically for large J. To find this behaviour 
we calculate /"(Aq) = —2/3. Also, 5'(Ao) = —2iJ/n. Thus, the leading steepest-descent 
calculation shows that ljj{2) diverges as J ^ oo: 

u;j(2) ~ / dXe'''''/' ~ (J ^ oo). (26) 



This reproduces the result of the Richardson extrapolation in ()14|) and we identify 
^ = 1.9544100476. . .. 

In summary, we have shown that as the number of energy levels increases, the 
normalised density of states /x(-E) approaches zero when E 1/2 and diverges when 
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E = 1/2. From this result and the normahsation condition satisfied by ^{E)^ we 
conclude that in this limit /i(-E) S{E — 1/2). Thus, according to the postulate 
that every quantum state associated with a given energy E must be realised with 
equal probability in microcanonical equilibrium, the density of states associated with 
a system having a nondegenerate linear energy spectrum approaches a delta function 
in the thermodynamic limit. It follows that in this limit the energy of the system can 
assume only one value. Whether an analogous result holds for an interacting system 
having a degenerate spectrum is an interesting open question. 
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